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Working in the hght-cone gauge, we find a simple procedure to calculate 
the autonomous one-loop Q^ evolution of the twist-three part of the nucleon 
gT{x,Q'^) structure function in the large-A^c limit. Our approach allows us 
C^ ', to investigate the possibility of a similar large-A'^c simplification for other 

^^ \ higher-twist evolutions. In particular, we show that it does not occur for the 

^ ' twist-four part of the fi{x,Q'^), g3{x,Q'^) and h^{x,Q'^) distributions. We 

^^ ■ also argue that the simplification of the twist-three evolution does not persist 

Q^ . beyond one loop. 
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Feyninan's parton model of incoherent parton scattering provides a transparent picture of 
what happens in a broad class of high-energy scattering processes. Modulo field theoretical 
logarithms, the parton model can be derived in quantum chromodynamics (QCD) in the 
form of factorization theorems [jl|. Better yet, QCD allows us to go beyond the naive parton 
model by consistently including the effects of the parton transverse momentum and coherent 
parton scattering. A simple example of coherent parton scattering is the interference of a 
single quark with a quark and a gluon in a nucleon target. To describe this phenomenon, it 
is necessary to introduce a three-parton light-cone correlation function 

M'^(a;,2/,Q2) = / ^^e'^^e^^(^-^)(P5|7A(0)zZ^"(/in)^(An)|P5) , (1) 

J Zn zn 

where n is a light-cone vector, tp a quark field, and \PS) the nucleon state. The general 
parton correlations involve more than one Feynman variable, and hence their scale (Q^) 
evolution is more complicated than the usual Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equations for the Feynman parton densities. Technically, the com- 
plication arises from the so-called higher-twist part of the correlations. 

Experimental study of parton correlations is challenging for a number of reasons. One is 
the lack of processes in which all Feynman variables in a parton correlation can be kinemat- 
ically controlled. For instance, in polarized lepton- nucleon deep-inelastic scattering (DIS), 
one can measure the structure function grix, Q^)- In the Bjorken limit, grix-, Q"^) is related 
to a y-moment of the above correlation function. Since a moment of Ma{x, y, Q^) does not 
evolve autonomously, knowing the entire gxiXyQ"^) at one scale is not sufficient to calculate 
it at another. This makes an analysis oi gxix^Q"^) data at different scales difficult. 

Several years ago, Ali, Braun, and Hiller (ABH) [§] made a remarkable discovery that 
in the limit of the large number of color Nf.., the twist-three part of gT{x,Q'^) does evolve 
autonomously at the one-loop level. The result has since been widely used in model cal- 
culations and analyses of experimental data |^. More recently similar results have been 
found for the evolutions of other twist-three functions /ii,(x, Q^) and e(x, Q"^) 0. Given the 
practical importance of the ABH result, a deeper understanding of the large Nc simplifica- 
tion is clearly desirable. Moreover, it is interesting to investigate the possibility of a similar 
simplification at two or more loops and for analogous twist-four correlations. 

In this paper we calculate directly the large- iVc evolution of grix, Q^) in the light-cone 
gauge. We find that the autonomy of the twist-three evolution arises from a special property 
of one particular Feynman diagram. Since this property is independent of the 7-matrix 
structure of the composite operators inserted, the ABH result generalizes immediately to 
the twist-three parts of hi{x,Q^) and e{x,Q^). Unfortunately, for various reasons we shall 
explain, there is no similar large- A''^ simplification for twist-four functions, nor for g2{x,Q'^) 
beyond one loop. 

We begin our discussion with a brief introduction to the grix, Q"^) structure function of 
the nucleon. In inclusive DIS, all information about the nucleon structure is summarized in 
the following hadron tensor, 

W^^(P,^,g) = ^ j d'ie^'^-^iPSWJ.iO.JMWPS) , (2) 

where J^ = J2q e^^gT^^g is the electromagnetic current and q is the spacelike virtual photon 



momentum. The antisymmetric part of the hadron tensor, M/t^'^l, is polarization-dependent 
and can be characterized in terms of the two structure functions gi{xB, Q^) and g2{xB, Q^)'- 

^[H ^ _,,,^a,^^ L^giix^ ^ ^^^^ _ ^^ q)p^fA^^\ ^ (3) 

where we have chosen the kinematic factors so that gi{xBiQ'^) and g2{xB,Q'^) survive the 
scahng hmit Q^ = —q^ — > cxd, u = P ■ q —^ oo and Xb = Q'^/'^u = finite. In Feynman's 
parton model, gi{xB, Q^) is related to the parton helicity density Aqa{x, Q"^) 
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:XB,Q') = lT.^l[^1ai^B,Q') + Aqa{-XB,Q')] , (4) 



where Ca is the electric charge and a sums over light quark species. 

The structure function g2{x b ■, Q"^) ■, however, does not have a simple parton model inter- 
pretation. Defining gT^XB^Q'^) = 5'i(a^_B,Q^) + 5'2(a^B, Q^), an operator-product-expansion 
analysis yields p 

^T(a:B,Q') = -5]e2(AgTa(xB,g') + AgTa(-a;B,g')) , (5) 

^ a 

where we have neglected all power and radiative corrections and 

AgTa(x,Q') = ^/^e^"^(P^x|^a(0)7^75^a(An)|PS^) . (6) 

The trouble with a parton model interpretation of AgT(j(x, Q^) can easily be seen in light- 
front quantization in which only the "good" component of the Dirac field ip^ = P^ip has a 
simple Fock expansion {P± = 7'^7^/2, 7^ = (7" ± 7^)/-\/2), whereas the "bad" component 
^_ = P^ip is constrained by the following equation of motion 

^_(An) = -^-^ ^tP^i\n)MXn) . (7) 

[In some sense ip- represents a quark-gluon composite.] Unlike Aqa{x,Q'^), AqTa{x,Q'^) 
contains a bad component because of the 7"*". 

For the same reason, the scale evolution of AqTa{x, Q^) is now more intricate than that 
of Aqa{x, Q"^). Its ra-th moment is written 

f^AqTaix,Q')x-dx = ^n,, ■ ■ ■n,„{PS^\e^^^--^-^\PS^), (8) 

where Q'^'y^^^^^^^^■^) = tjj'j'^iD^^^'^ ■ ■ ■iD^^"^ip, with [jii-'-jjin) indicating symmetrization of the 
indices and removal of the traces. The 6'-operator contains both twist-two ^(°"^i ■■■^") (totally 
symmetric and traceless) and twist-three ^[°"(/^i]^2-Mn) (mixed symmetric and traceless) con- 
tributions, where [cr/ii] denotes antisymmetrization. It turns out, however, that for a given 
symmetry structure there are multiple twist-three operators. In fact, a complete basis of 
these operators was first identified in H, 



S^ = ^piD^f"^ ■ ■ ■ iDf'^-'gF''''HD''^+' ■ ■ ■ iDf'^-^Y^^ip , (9) 

where i = 1, ..., n — 1. The operator ^['^(^il^a ■A'n) jg j^g^ g^ special hnear combination of them, 
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^['^('^^i^-'^") = ^^;^^ E(^ - ^)(^r - K-^ + sr + s:_,) . (lo) 

The anomalous dimension matrix in the above operator basis was first worked out by 
Bukhvostov et al. and later reproduced by a number of authors with different methods M . 
The result is what one would generally expect. To evolve the matrix element of ^['^(/^il'^2--Mn)^ 
it is not enough just to know it at an initial scale — one must know all the matrix elements 
of W^^ = R^- Rl_i + ^r + Sl_i there. 

By studying the anomalous dimension matrix in the large N^ limit, Ali, Braun and 
Hiller found that the eigenvector corresponding to the lowest eigenvalue is just the linear 
combination of twist-three operators on the right-hand side of Eq. (0). In other words, the 
twist-three part of Aqxaix^Q"^) evolves autonomously in this limit. To better understand 
ABH's result, we calculate the large- iVc evolution of AgTa(a^, Q^) directly. We start with the 
mixed-twist operator ^°"(MiA'2-Mn) j^ gq_ (^g^ g^j^j look for possible divergences when inserted 
in multi-point Green's functions. To reduce the number of Feynman diagrams, we choose 
the light-cone gauge A^ = and take the ± + ■ ■ ■ + component of the ^-operator. Let's call 
the resulting operator 6'„ = ilj'j-^'j5{id'^)""ip, and its twist-two and twist-three parts 9n2 and 
9n3, respectively. The Feyman rule for 9n is simply 7"'"7^(fc"'")", where k is the momentum of 
the quark. 

By light-cone power counting, we need only consider two- and three-point functions. 
Since the external lines carry color, we must ask what type of diagrams dominates the 
large Nc limit. The simple rule we find is that when all external lines are drawn to one 
point (infinity), the planer diagrams are leading. All one-particle-irreducible (IPI) leading 
diagrams with one 6 insertion are shown in Fig. 1. 

The ultraviolet divergences in the two point Green's function can obviously be subtracted 
with the matrix element of 6n itself. The only diagram in which the divergences may not be 
subtracted by On is Fig. lb. An explicit calculation shows that the ultraviolet divergences 
correspond to the following local operator: 






-CA^lnQ 



1 n— 1 



(^ + 2)^0 

(11) 



/n+l 1 1 \ _ 

+ ( E - - ^(^^^ ) (^^^ ^7^75(^5+)"- V + V^(^9+)"-S^75 l^^P±i^ 

where we have neglected the contributions of light-cone singularities which will be cancelled 
eventually. Notice that the first term is present in the twist-two operator 

en2 = -^ fv'7^75(^5+)> + E ^7^75{^^^y^D^{^^+r-'-'^^^ , (12) 

^ + -•- V i=0 / 
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(a) (b) 

FIG. 1. Two and three-point IPI Feynman diagrams contributing to the evolution of On in the 

large Nc limit. 

and the remaining two terms can be converted to On by using the equation of motion in 
Eq. (0). Thus we easily arrive at the ABH conclusion that ^„3 evolves autonomously in the 
large- A'"c limit. 

Including the contribution from Fig. la as well as the one-particle-reducible ones that 
cannot be neglected in the light-cone gauge, we obtain the following equation, 



ddn as{Q'')N, 



rflnQ2 27r 



n + 2 \ p' 2 




(13) 



Separating out the twist-two and twist-three parts, we not only recover the well-known 
twist-two evolution, but also the twist-three result 

i?^.-m(_2i:'i, 1 a),„3. (14, 

hiQ-^ Ztt \ ^2 n + 1 2 J 

which is identical to the result in Ref. 0. 

It is quite clear that the i-independence of the coefficients in the sum of Eq. ( pT]) is the 
key for the autonomous evolution of 6n3- On the other hand, this property is not totally 
unexpected if one inspects Fig. lb more closely. Interpreting this diagram in the coordinate 
space, we see that the internal gluon propagates homogeneously from one quark to the 
other. By homogeneously, we mean that at any point along the path of the propagation, 
the gluon behaves exactly the same way, except, of course, at the points where the gluon 
and quarks interact. Now the spatial location of the interaction with the external gluon 
determines the number of derivatives before and after the gluon field in the subtraction 
operator. Since the internal gluon propagation is homogenous, different locations of the 
triple-gluon vertex should produce similar physical effects. Therefore, the coefficients of the 
subtraction operators tjj fi'j5{id~^yiD-^{id~^)'^~^~^tp should be independent of z. On the other 
hand, the two extra terms in Eq.(pd]) correspond to the triple-gluon vertex just next to the 
external quark lines, where the homogeneity is lost. 

Since the homogeneous property of the internal gluon line is independent of the gamma 
matrix structure of the operator inserted, we conclude that the other twist-three distributions 
e{x,Q'^) and hi^x^Q"^) evolve also autonomously in the large Nc limit. A quick calculation 
confirms the evolution equations found in Ref. @]. 



Encouraged by the success of the above approach, we apply it to the analogous twist-four 
evolution. In Ref. ^, the three one-variable distributions f^iXjQ'^), g^i^x^Q"^) and h^i^x^Q"^) 
are shown to contain twist-four. For example, fi{x) is defined as 

U^) = ^ J ^{P\mi-H^n)\P) . (15) 

It was shown in Ref. that f4{x,Q^) contributes to the l/Q^ term of the longitudinal 
scaling function F^ of the nucleon 

Fl{xb, Q') = ^^ E ea/4a(xB, Q') , (16) 



a 



where we have neglected higher-order radiative corrections. Here, autonomous evolution of 
f4:{x, Q"^) would simplify the analysis of Fl data immensely. 
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FIG. 2. Four-point IPI Feynman diagrams contributing to the evolution of O in the large N^ 
limit. 

In the large N^ limit, we consider one insertion of the operator O = ip^^ {id^)^ip into 
two-, three- and four-point Green's functions. At one-loop order, the IPI two- and three- 
point diagrams are identical to those in Fig. 1 and the IPI four-point diagrams are shown 
in Fig. 2. Only the three and four point diagrams can potentially destroy the autonomous 
evolution of O. Let us start with Fig. 2a. One of the divergent contributions from this 
diagram introduces the following local subtraction 

E i^iV>± fiiid^yi]/)i_ {id+Y~"^ + h.c. (17) 

i 

where all the coefficients are independent of i again because of the homogeneity of the gluon 
propagator. Using the equation of motion, we can write this as 

Y^i^iid+Yi p^iid^y'-'-^ij + h.c. (18) 

i 

Since this operator cannot be reduced to either the twist-two or twist-four part of O, the 
evolution of the latter cannot be autonomous unless this contribution is cancelled by other 
diagrams. The only other diagram containing the same divergence structure is Fig. lb with 
an insertion of O. Unfortunately, our explicit calculation did not produce this cancellation. 
The same phenomenon occurs for the twist-four part of gs^x, Q^) and /i3(x, Q^). 



Thus, the large Nc simphfication seems to happen only for the evolution of the twist-three 
part oi gxix, Q'^), hii^x, Q^) and e(x, Q^). Does it happen for them at two and higher loops? 
In Fig. 3, we show two examples of Feynman diagrams that we suspect break the autonomy 
of the ^3„-evolution, i.e., they may contain divergences that cannot be subtracted by 9n2 and 
6n3 only. Our suspicion is based on the inhomogeneity of the gluon progator. The internal 
gluon that propagates from one quark to another has different wavelengths in the different 
parts of the propagation. Its interaction with the external gluon is different at different 
spatial locations. Thus the subtraction operators have different coefficients depending on 
the number of derivatives before and after the external gluon field. An explicit calculation 
of Fig. 3a confirms our suspicion. 





(a) (b) 

FIG. 3. Some two-loop IPI Feynman diagrams that might break the autonomy of the On in the 

large N^ limit. 

This leaves us with only one possibility for autonomous two-loop evolution of 6n3'. the 
unwanted structures cancel in the sum of all large- iVc two-loop diagrams. Calculating all 
those diagrams is a big task. However, even without an explicit calculation, we do not 
expect the cancellation to happen. The fundamental reason is that large A*";, represents only 
a selection of a subset of Feynman diagrams, whereas the result of an individual diagram 
is independent of the large- iVc limit. Cancellations of a structure do not happen among 
Feynman diagrams unless there is a symmetry. 

Therefore we conclude that the autonomy of one-loop evolution for a set of special twist- 
three distributions at large Nc seems accidental. In the light-cone gauge, it can be easily 
traced to a special property of Fig. lb. The simplification does not happen for the analogous 
twist-four distributions at one loop, nor for those twist-three distributions at two or higher 
loops. Nonetheless, the discovery of Ali, Braun, and Hiller remains as a significant step 
forward in the study of the g2{x, Q^) structure function. Without the autonomous one- loop 
evolution, an analysis of experimental data on the twist-three contribution would be severely 
constrained. 
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